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Interaction between Faraday rotation and Cotton- Mouton effects in 
polarimetry modeling for NSTX 
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The evolution of electromagnetic wave polarization is modeled for propagation in the major radial direction 
in the National Spherical Torus Experiment (NSTX) with retrorcflcction from the center stack of the vac- 
uum vessel. This modeling illustrates that the Cotton-Mouton effect-elliptization due to the magnetic field 
perpendicular to the propagation direction-is shown to be strongly weighted to the high-field region of the 
plasma. An interaction between the Faraday rotation and Cotton-Mouton effects is also clearly identified. 
Elliptization occurs when the wave polarization direction is neither parallel nor perpendicular to the local 
transverse magnetic field. Since Faraday rotation modifies the polarization direction during propagation, it 
must also affect the resultant elliptization. The Cotton-Mouton effect also intrinsically results in rotation of 
the polarization direction, but this effect is less significant in the plasma conditions modeled. The interaction 
increases at longer wavelength, and complicates interpretation of polarimetry measurements. 
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I. INTRODUCTION 

Polarimetry is a powerful technique for probing mag- 
netic field equilibria and fluctuations, plasma density 
and cur ren t density profiles in magnetically confined 
plasmasP^ It measures changes in the electromagnetic 
(EM) wave polarization caused by propagation through 
a magnetized plasma. It is routinely used on conventional 
high aspect ratio tokamaks (e.g. JET^l) and reversed field 
pinches (REP, e.g. MST 4 ). However, no detailed study 
of polarimetry has been performed for propagation in 
the major radial direction in spherical tori. In contrast 
with conventional tokamaks and RFPs, in spherical tori 
both magnetic field strength and direction vary strongly 
in the major radial direction. The modeling is moti- 
vated to guide the design of a polarimeter system planned 
for the National Spherical Torus Experiment (NSTX)P. 
It calculates the evolution of EM wave polarization in 
the major radial direction with retroreflection from the 
center stack of the vacuum vessel. The Cotton-Mouton 
effect-elliptization due to the magnetic field perpendicu- 
lar to the propagation direction-is shown to be strongly 
weighted to the high-field region in NSTX. An interac- 
tion between the Faraday rotation and Cotton-Mouton 
effects is also clearly identified. Elliptization occurs when 
the wave polarization direction is neither parallel nor per- 
pendicular to the local transverse magnetic field. Since 
Faraday rotation modifies the polarization direction dur- 
ing propagation, it must also affect the resultant elliptiza- 
tion. The Cotton-Mouton effect also intrinsically results 
in rotation of the polarization direction, but this effect 
is less significant in the plasma conditions modeled. The 
interaction is shown to increase with wavelength. 

The interaction is present when the magnetic field has 
both parallel and perpendicular components with respect 
to the wave propagation direction. It complicates the in- 
terpretation of polarimetry measurements, especially at 
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longer wavelength. Previous polarimetry studies focused 
separately on Faraday rotation^ or the Cotton-Mouton 
effeclP. Recent results combining measurement and mod- 
eling including both effe cts on JET also assume one effect 
or the other is smallJ^^ However, this will not be the case 
in ITER, where both effects are large.^ 

In the following sections, the polarimetry model is 
described and results are shown for modeling using a 
plasma density profile from Thomson scattering measure- 
ment and a magnetic field profile from EFIT in NSTX. 



II. POLARIMETRY MODELING DESCRIPTION 

Several assumptions are made to simplify the calcula- 
tion of polarization evolution. A cold plasma model is 
adopted, which excludes corrections from finite temper- 
ature effects^. The plasma is assumed to be collision- 
less, so the beam experiences no dissipation. The WKB 
approximation^ is used, i.e. plasma parameters are as- 
sumed to be slowly varying (\B\ ^> \(l/k)(dB/dz)\,n ^> 
|(1/ k)(dn/dz)\), so cutoffs and resonances are not consid- 
ered. Only the electron response is included; the contri- 
bution from ion motion is ignored (oj p i,u) C i <C uj pe ,uj ce <C 
ui, where w is the circular frequency of the probing beam, 
uj p = y / nq 2 /meo, w c = \q\B/m, subscripts '£' and 'e' 
stand for ions and electrons respectively). Refraction 
is also neglected, so it is assumed that the beam path 
through the plasma is straight and that the beam nei- 
ther diverges nor converges. 

The polarization of an EM wave evolves as it propa- 
gates through a magnetized plasma due to plasma bire- 
fringence. The plasma features a pair of fast and slow 
characteristic modes (i.e. EM waves that propagate with 
their polarizations unchanged) whose phase velocities 
and polarizations are determined by local plasma param- 
eters at any position along the wave path [Fig. Ilia)]. 
The two modes are generally elliptically polarized with 
orthogonal polarization directions and opposite handed- 
ness. The fast mode has a major axis perpendicular 



fast mode. 




FIG. 1. (a) Fast (solid line) and slow (dashed line) character- 
istic modes. The fast mode has a major axis perpendicular to 
B± and right-handedness with respect to By. '_!_' and '||' are 
defined with respect to the propagation direction, (b) Polar- 
ization properties are characterized by \ [elliptization angle, 
where right/left handedness are represented by +/— sign) and 
ip [polarization direction angle); their ranges are also shown. 



to B± and right-handedness with respect to By. ('_L' 
and '||' are defined with respect to the propagation di- 
rection) An EM wave of any polarization may be rep- 
resented as some combination of this pair of character- 
istic modes. The combined polarization is sensitive to 
the relative phase of its two components, so a difference 
in their phase velocities causes the polarization to evolve. 
The Faraday rotation and Cotton- Mouton effects are two 
well-known special cases. In Faraday rotation, where the 
wave propagates parallel to a magnetic field, the charac- 
teristic modes are circularly polarized. For the Cotton- 
Mouton effect, where the wave propagates perpendicular 
to a magnetic field, the fast and slow modes, which are 
linearly polarized, are extraordinary and ordinary modes 
respectively. 

The modeling uses the Mueller-Stokes calculus^. The 
ellipse representing the polarization of a single-frequency 
EM wave is characterized by two parameters, cllipti- 
zation angle, x^ an d polarization direction angle, tp. 
[Fig. [lib)] The polarization state of the wave can be 
mapped by the Stokes vector s (Eq. nl to a point on 
a unit sphere in an abstract space referred to as the 
Poincare sphere. For instance, alignment of s with the 
S3 axis corresponds to circular polarization, while a van- 
ishing S3 component corresponds to linear polarization in 
the laboratory frame. As a wave propagates through a 
magnetized plasma, its polarization evolves and the cor- 
responding Stokes vector traces out a trajectory on the 
Poincare sphere. Each small step of the trajectory results 
from a small rotation of the Stokes vector around an axis 
given by the vector fi which is determined everywhere 
along the wave path by local plasma parameters. (Eq. [2]) 
The z coordinate indicates position along the wave path 
and c is speed of light in vacuum. 



ds(z) 

dz 



Q(z)xi?(z),Q 



, , 2 1 , 2 



2cw{uj 2 - W 2 J 



{Bl-BD/B 
2B x B y /B~ 

2(; 



:)B Z /B 

(2) 
The origin of the interaction between the Faraday ro- 
tation and Cotton-Mouton effects can be seen clearly 
from the preceding geometrical description of the polar- 
ization evolution. Both the Faraday rotation and Cotton- 
Mouton effects are directly related to the components of 
ft. A non-vanishing O3 gives rise to Faraday rotation by 
causing a rotation of s about the S3 axis. This corre- 
sponds to a rotation of the wave polarization ellipse in 
the laboratory frame. A non-vanishing f2i or SI2 gives 
rise to the Cotton-Mouton effect by causing a change in 
S3 and therefore in the ellipticity of the wave polarization 
ellipse. However, the way in which S3 changes clearly de- 
pends on the direction of s relative to f2, which can be 
influenced by Faraday rotation. 

The interaction may also be seen from the differential 
expressions relating changes in x an d ^ to the Faraday 
rotation (FR) and Cotton-Mouton (CM) effects^ 



d X = 2 Sin2 ^ d ^lcM 



(3a) 



d^ = dV>|FR - - tan 2x cos 2i))d5\cM (3b) 
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6 is the relative phase between the x and y components 
of the wave electric field. Eq. [3] assumes the coordinate 
system illustrat ed i n Fig. nib) where the x axis is aligned 



with B±. Eq. 3a shows the sensitivity of elliptization 
to ip. Faraday rotation modifies ip, so it affects elliptiza- 
tion. Eq.|3b|shows that the Cotton-Mouton effect also in- 
trinsically causes polarization rotation, but for the cases 
modeled here, this contribution to the total polarization 
rotation proves to be relatively small. This is because 
through out the majority of the wave path, either the 
absolute elliptization angle is small (i.e. |x| <§C 45°) or 
the Cotton-Mouton effect is overwhelmed by Faraday ro- 
tation (i.e. |(W|cm| "C |c1-0|fr|)- The following discussion 
will focus on the impact of Faraday rotation on elliptiza- 
tion. 



III. MODELING RESULTS AND DISCUSSION 
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Many of the modeling results discussed here are ob- 

(1) tained for a 288 GHz (A = 1.04 mm) probing microwave 

in a base case plasma (shot # 124764, 0.325 sec) that 
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FIG. 3. Elliptization angle (x) evolution along chord in mid- 
plane for waves with horizontal linear polarization at launch 
(i.e. in toroidal direction). Vertical solid line indicates plasma 
center (i.e. peak density). The mirror is mounted on the cen- 
ter stack. (/ = 288 GHz, X = 1.04 mm) 



FIG. 2. (a) Plasma density profile of base case for modeling, 
(shot # 124764, 0.325 sec), (b) Toroidal (solid line) and ver- 
tical (dashed line) magnetic field along major radius in the 
mid- plane (they vary little with height near the mid-plane). 
(c) Horizontal (i.e. radial) magnetic field along major radius 
0.1m above (blue dashed line), below (green dashed line) and 
in the mid- plane (red solid line). 



represents conditions likely to be encountered in NSTX. 
The base case is an L-mode plasma [Fig. [2] with a major 
radius of R = 0.85 m and a minor radius a = 0.67 to. 
The density profile is centrally peaked with a maximum 
of n = 4.7 x 10 19 to~ 3 at R axis = 1.0 to. The elec- 
tron plasma and cyclotron frequencies are 61.4 GHz and 
10.5 GHz on axis. The millimeter wavelength used is 
longer than typical for polarimetry systems, but it is 
a good compromise between two competing constrains 
for the base case. At longer wavelength the effects of 
the plasma on polarization are stronger, allowing for 
more sensitive measurement of magnetic fluctuations. At 
shorter wavelength, refraction becomes less significant. 

The modeling shows that the magnitude of the ellipti- 
zation angle increases most rapidly when the wave is in 
the high- field region (R < R ax is)- [Fig. |3] This stands 
in contrast with conventional tokamaks. Elliptization is 
sensitive to the strength of the perpendicular magnetic 
field (Eq. 3a), of which the toroidal magnetic field Bt 
is a significant component. Bj> varies approximately in- 
versely with major radius in both conventional tokamaks 
and spherical tori, but in spherical tori, the variation is 
much stronger because of their relatively low aspect ra- 
tio. For instance, in NSTX (Ro/a ~ 1.27) Bt varies from 
0.2 T at the outer edge (R = 1.6 m) to 2 T close to the 
center stack (R = 0.2 to). [Fig. (2jc)] 

Modeling shows that the evolution of the elliptization 
depends strongly on the polarization direction in the 
high-field region. Fig. [4] shows the dramatically differ- 
ent elliptization evolution for two waves launched in the 
mid-plane with launch angles of 0° and 45° . This depen- 
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FIG. 4. Elliptization angle (x) evolution for horizontal (solid 
line) and 45° (dashed line) linear polarization at launch along 
chord in mid-plane. Zero elliptization (i.e. linearly polarized) 
is highlighted on the grid. (/ = 288 GHz, X = 1.04 mm) 



dence is expected since elliptization is strongly weighted 
to the high- field region and sensitive to ip (Eq. 3a). For a 
chord in the mid-plane, the polarization direction in the 
high-field region is determined by the launch angle since 
Faraday rotation is very weak there (By is weak in the 
mid-plane) . 

Of particular interest, the modeling shows that Fara- 
day rotation can play a significant role in elliptization. 
Chords away from the mid-plane can have significant Bu , 
so Faraday rotation can substantially change the polar- 
ization direction of the wave before it enters the high- 
field region. Fig. [5] compares the elliptization evolution 
of a wave launched with horizontal linear polarization 
both with and without the influence of Faraday rotation. 
The final elliptization of the wave is very different for the 
two cases. The modeled chord is 0.1 to above the mid- 
plane, where |Bm| reaches a maximum of 0.024 T along 

the chord. For the case without Faraday rotation, Bu is 
simply set uniformly to zero along the chord. The im- 
pact of Faraday rotation on elliptization identified here 
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not only in the plasma modeled here, but whenever the 
wavelength is sufficiently long regarding the plasma con- 
ditions, such as in the planned 118 fmi polarimeter in 
ITER 11 . Under these conditions, the approximate ex- 
pressions (integral forms of Eq. El for the polarization 
rotation and elliptization are not appropriate. Also, the 
interpretation of an array of chord measurements used to 
characterize the equilibriu m 16 ! 1 I becomes more compli- 
cated. The profile of final polarization direction versus 
chord impact parameter is affected by the interaction, 
leading to a change in both the zero crossing and slope. 



FIG. 5. Elliptization angle (x) evolution of wave with hori- 
zontal linear polarization at launch along a chord 0.1m above 
mid-plane with (solid line) and without(dashed line) Faraday 
rotation. Faraday rotation is eliminated by setting Bu — 
along chord). (/ = 288 GHz, X = 1.04 mm) 
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FIG. 6. Difference between the final elliptizations with and 
without Faraday rotation (Ax) versus wavelength for different 
launch angles. Ax is normalized by the difference between the 
maximum and minimum final values of x without Faraday 
rotation. 



IV. SUMMARY AND CONCLUSION 

This work models the evolution of EM wave polariza- 
tion along major radial chords, with retroreflcction on 
NSTX, using Mueller-Stokes calculus. Most of the model- 
ing focuses on 288 GHz microwaves launched with linear 
initial polarization in a base case plasma that represents 
conditions commonly encountered in NSTX. The mod- 
eling shows that the Cotton-Mouton effect is strongly 
weighted to the high-field region of NSTX. An interac- 
tion between the Faraday rotation and Cotton-Mouton 
effects is also clearly identified. Elliptization occurs when 
the wave polarization direction is neither parallel nor per- 
pendicular to the local transverse magnetic field. Since 
Faraday rotation modifies the polarization direction dur- 
ing propagation, it must also affect the resultant ellipti- 
zation. The Cotton-Mouton effect also intrinsically re- 
sults in rotation of the polarization direction, but this 
is less significant to the modeling results presented here. 
The interaction identified here is shown to increase in 
significance with wavelength. Care has to be taken in 
interpreting the polarimetry measurement if both effects 
are large. 



is a primary element of the interaction between the two 
effects. 

Modeling shows a significant wavelength dependence in 
the strength of the impact of Faraday rotation on ellip- 
tization. Both the Faraday rotation and Cotton-Mouton 
effects are expected to become stronger with increasing 
wavelength (i.e. lower frequency). (Eq. [4]) However, it is 
not obvious whether the impact of Faraday rotation on 
elliptization should become more or less significant as the 
wavelength increases. To assess this, the change in the 
final value of \ caused by including Faraday rotation is 
calculated. [Fig. [6] The change A\ is normalized by the 
difference between the maximum and minimum final val- 
ues of x without Faraday rotation. This normalization 
factor serves as a measure of the strength of the elliptiza- 
tion effect. Fig.|6]shows that the relative effect of Faraday 
rotation on elliptization increases with wavelength. 

The interaction complicates the interpretation of po- 
larimetry measurements if both the Faraday rotation 
and Cotton-Mouton effects are large. This is an issue 
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